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^^ i Abstract. In this paper we are concerned with the global existence of smooth 

^~> , solutions to the turbulent flow equations for compressible flows in R 3 . The 

Q\ ■ global well-posedness is proved under the condition that the initial data are 

close to the standard equilibrium state in i? 3 -framcwork. The proof relies 

on energy estimates about velocity, temperature, turbulent kinetic energy and 
P^ , rate of viscous dissipation. We use several new techniques to overcome the 

difficulties from the product of two functions and higher order norms. This is 

the first result concerning k-e model equations. 
^ ' AMS Subject Classification 2000: 35Q35, 35A01, 76F02. 
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/"5 ' 1. Introduction 

in 

l/~) . We consider in this work the turbulent flow equations for compressible flows on 

en 

o 

(N 



l.i 



' p t +div(pu) = 0, 

(pu)t + div(pu ® u) — Au — Vdivu + Vp = — |V(/9fe), 

(ph) t + div(puh)-Ah=^+S k , 

{pk) t + div(puk) — Afc = G — pe, 

(pe) t + div(pue) - Ae = ^ - ^pl, 
^ (p, u, h, k, e)(x, *)|t=o = (Po(x),u (x),ho(x),ko(x),E Q (x)), 

With S k = [m(J^ + ¥-)- &ijli%£]Tr- + ^^7^' G = ^[Peilr- + W L )~ 
K "\Q X j oxi ' 3 L J " oxh J oxj p z axj oxj 1 oxj^^^oxj oxi ' 

\bij{pk + /i e |^)], where % = if i ^ j, 5 V] = 1 if i = j, and p, p t , Ci, C 2 are 
four positive constants satisfying p + p t = Pe- 

Here p, w, h, k and e denote the density, velocity, total enthalpy, turbulent 
kinetic energy and rate of viscous dissipation, respectively. The pressure p is a 
smooth function of p. In this paper, without loss of generality, we have renormalizcd 
some constants to be 1. 
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The system (jl.lj) is formed by combining effect of turbulence on time-averaged 
Navier-Stokes equations with the k-e model equations. Turbulence stands out as a 
prototype of multi-scale phenomenon that occurs in nature. It involves wide ranges 
of spatial and temporal scales which makes it very difficult to study analytically 
and prohibitively expensive to simulate computationally. Up to now, there are not 
any general theory suitable for the turbulent flows. Many, if not most, flows of 
engineering significance are turbulent, so the turbulent flow regime is not just of 
theoretical interest. Hence, research about the above system (|1.1[) is initial and 
very important. 

There are many results about the simple model, such as classical Navier-Stokes 
equations and incompressible Navier-Stokes equations with variable density. 

For classical Navier-Stokes equations (let p = C, k = e = 0), there are many 
authors doing research about them, but the uniqueness and regularity of weak 
solutions for general initial data arc open. In particular, Fujita and Kato [1 61 117] 
proved the local well-posedness of 3D Navier-Stokes equations for large initial data 
and the global well-posedness for small initial data in the homogeneous Sobolcv 
space H?. If taking uq small with respect to the viscosity in the Lebesgue space 
L N (M. N ), the classical Navier-Stokes equations are globally well-posed (see p"8ll22j ). 
This result has been adapted to the case of bounded domains by Y. Giga and T. 
Miyakawa in [19] and to the case of exterior domains by Y. Giga and H. Sohr in 
[20] . and H. Iwashita in [21]. In the half-space case, the well-posedness issue has 
been studied by H. Kozono in [53] (see also [B] for results related to critical Bcsov 
spaces with negative index of regularity). In the absence of heat conduction, it 
was proved by Z. Xin that any non-zero smooth solution with initial compactly 
supported density would blow up in finite time(see |30j). As a reasonable starting 
point, we will therefore restrict our work to solutions such that p remains positive. 

For incompressible Navier-Stokes equations with variable density (let k = e = 0), 
there are also many results. Specially, O. Ladyzhcnskaya and V. Solonnikov [53] 
studied the existence of strong (smooth) solutions with positive density whereas 
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the theory of global weak solutions with finite energy has been performed in the 
book |^5]. Abidi Q] and Danchin-Mucha [HJ- [XO] have obtained several results about 
incompressible Navicr-Stokes equations with variable density in the critical spaces. 
For more about well-posedness in other critical spaces, readers can refer to [21 [3j [4j 

EnzumnaiiainiiiaEBi. 

For the system (jl.ljl . there is not any result. This paper is devoted to the study 
of the global existence of smooth solutions for the system (jl.ll) under suitable 
assumptions. Our result is expressed in the following. 

Theorem 1.1. Assume that initial data are close enough to the constant state 
(p, 0, 0,fc, 0), i.e. there exists a constant 6q such that if 

|| (p - p,u ,h ,k a - fc,£o)||ff3( R 3) < Sq, (1.2) 

then the system il.l]) admits a unique smooth solution (p, u, h, k, e) such that for 

any t S [0, oo), 

\\(p-p,u,h,k-k,e)\\ 2 H3 + / \\Vp\\ 2 H 2 + \\(Vu,Vh,Vk,Ve)\\ 2 H3 ds 



'0 

where C is a positive constant. 



< C\\(p - p,u ,h ,k - fc,£o)Hi 



Remark 1.2. The existence of the local solution for (|1.1|) can be obtained from the 
standard method based on the Banach theorem and contractivity of the operator 
defined by the linearization of the problem on a small time interval (see also [27j [28j 
[29]). Hence, we omit the local existence part for simplicity. The global existence 
of smooth solutions will be proved by extending the local solutions with respect to 
time based on a-priori global estimates. 

Remark 1.3. We use several new techniques to overcome the difficulties from the 
product of two functions and higher order norms. In a separate paper we consider 
the convergence rate of smooth solution to the constant state. 

Notation: Throughout the paper, C stands for a general constant, and may 
change from line to line. The norm || (A, B)\\x is equivalent to \\A\\x + \\B\\x- The 
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notation L P (R 3 ), 1 < p < oo, stands for the usual Lebesgue spaces on M 3 and || • || p 
denotes its LP norm. 

2. A PRIORI ESTIMATES 

In this section we establish a-priori estimates of the solutions, we assume that 
(v, u, 6, k, e) is a smooth solution to (|1.1|) on the time interval (0,T) with T > 0. 
We shall establish the following proposition. 

Proposition 2.1. There exists a constant <5 <^C 1 such that if 



sup \\(p-p,u,h, k-k, e)\\ H 3 < 5, 

0<t<T 

then for any t € [0,T], there exists a constant C\ > 1 such that it holds 
\\(p-p,u,h,k-k,e)\\ 2 H3 + f ||Vp||^ 2 + ||(Vu,V/i,Vfc,Ve)||^3rfs 



(2.1) 



(2.2) 



< Ci||(Po - p,uo,ho,k -fe,£ )||lf3. 



Proof. First, let p = a + p, k = m + k, f'(p) = S-iEl ) we rewrite the system (|1.1[) as 



follows: 



at + div((a + p)u) = 0, 

u t + u-Vu- -^{Au + Vdivu) + V[/(a + p) - f(p)] 

= -3(iFp)V((a + p)(m + k)), 

ht + u-Vh- ^Afc = -f\a + p)(a + p)divu + ^S k , 

mt + u • Vm - ^Am = ^G - e, 

s t + u ■ Ve - -4=Ae = , , c ^ Ge , ,, - ^fr, 

1 a+p (a+p)(m+k) m+k ' 



(2.3) 



t (a, u, a, 77i,£)(x,i)| t= o = (ao(x),«o(a:),/io(a:),mo(a;),eo(x)), 
with S k = [/i(fHl + |Hi) _ |^i^]|i + j^^-^r, G = fJMfj + f£) 



f<J«((o + /5)(m + A) + Me^r)]. 



From a priori assumption (|2. 1 [) and the Sobolev inequality together with the first 
equation of (|2.3[) . we have 



sup | (a, a*, Va, u, Vu, a, S7h,m, Vm, e, Ve)| < C||(a, u, /i, m, e)||#3 < C<5. (2.4) 



Moreover, 



-<p = a + p<2p, -<k = m + k<2k, 



(2.5) 



GLOBAL SOLUTION FOR TURBULENT FLOW EQUATIONS 5 

and 

< — < /'(/?) < C < oo, \f (n) (p)\<C for any positive integer n, (2.6) 
Co 

with Co a positive constant. 

In what follows, we will always use the smallness assumption of S and (|2.4[) - (|2.6I) . 
We divide the a priori estimates into three steps. 

Step 1: L 2 -norms of u, h, m, £ 

Multiplying the second equation of (|2.3p by u and integrating over R 3 , one can 
deduce that 

54|Nl2 + -^(||V«||l + ||divtt|||) 
l at a + p 

u ■ Vit • udx + / (/(a + p) — f (p))dwudx 

1 Jr3 (2.7) 
2 [(Va ® w) : Vu + divu(Va • u)]dx 



(o + p) 

2 

— — — V((a + p)(m + k)) ■ udx. 
j 3 3(a + p) 

With the help of Holder's inequality, VF 1 ' 2 (IR 3 ) <->• L 6 (R 3 ) and fl2l}, we estimate 

the right-hand side of ()2.7|) as 



u • Vu • udx < C||w|| 6 ||Vm|| 2 ||m|| 3 < C<5||Vu|||, (2.8) 

(2.9) 

V((o + p)(m + fc)) -lidx 
3(a + p) 

\7 n . ii — In J- n\\7 <m . i/Ht \ * ' 



/ - — - [( Va ® u) : Vu + divu( Va • u)]d 

J R3 [a + py 

< C*||Va||2||Vu|| 2 ||u|| L =o < C(5||(Va,Vu)j|2, 
2 



f 2(m + k) 2 

J R3 3{a + p) 3(a + p) 

< Ci5||(Va,Vm)|||, 



(/( a + P) - /(/o))divuda; 

= -/[/(a + p-)-/(p)](^^)d, 

Jr3 a + P 



. f F{a)dx - [ f '(P + e _ a K u-Vadx 
Jr3 Jr3 a + p 



(2.11) 



<~ / F(a)dx + C<5||Va||i, 
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with 9 £ (0, 1) and F(a) is denned as 

F(a)=r« a + d- f W d s. (2.12) 

Jo s + P 

Combining ([2~g]) - ([2~TTj) with (|27f)l . we get 

o-fA\Hl+ / i ? (a)da;] + — — (||V U ||i + ||diHll) < C(J||(Vo,V«,Vm)||§. (2.13) 

Multiplying the energy equation, governing equation for turbulent kinetic energy 
k and e-equation of (|2.3p by /i, to and e respectively, and integrating them over the 
whole space R 3 , similarly we can get that 

— 11*112 + r^llVfclli 

2 at a + p 

= / rr^-Va • V/i • /idx - / /'(a + p)(a + p)divu/ida; (2.14) 

J M 3 (a + pY J r3 

+ / -Sk ■ hdx — I u ■ V/i ■ /idx, 

J R 3 a + p J R 3 

vjl\\m\\Z + ——\\Vm\\l 

2 at a + p 

1 f f 1 

— Vm • Va • mda; — / e ■ mdx + / :G • mdx (2.15) 



( a + P) 2 Jm? J R 3 a + p 



u 


■ Vto 


■ mdx, 










2dr 


41 + 


a + p 


|Ve| 


2 
2 



Ve-Va-eda;+ / ^- ^dx (2.16) 



(a + p) 2 7 K 3 (a + /9)(m + fc) 

2 



CV 3 



dx — / it • Ve • £d:c. 



J3 TO + fc 

A direct computation gives that 

/ t — — ^Va-Vh-h<te:<C\\Va\\2\\Vhh\\h\\L<>° <C6\\(Va,Vh)\\%, (2.17) 
J R 3 (a + pY 

/'(o + p)(o + p)div«ftda;<<7||Vo|| 2 ||Vtt|| 2 ||ft||3<C<y||(Vo,V«)||l, (2.18) 

-Su ■ hdx 

a + p 

1 r, r^£_ + ^_\_^.. du k du l fi t dp da ^ (2.19) 

a + p dxj dxi 3 %3 dxk dxj (a + p) 2 dxj dxj 

< C5\\(Va,Vu)f 2 , 



u-Vh- hdx < C||w|| 3 ||V/i|| 2 ||/i|| 6 < Co||V/i|| 2 , (2.20) 
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/ - =^Vm- Va-mdx< a||Vm||2||Va|| 2 ||m||z,=o <C'5||(Vo,Vm)||| ] (2.21) 

J R3 (a + py 

e ■ mdx < C||e||2||m|| 6 ||m + fc|| 6 ||m + k\\ 6 < C5\\Vm\\j, (2.22) 



1 -_G- m dx 



a + p 



7Tn7^ [fIe{ 7^ + 7^ ] ~ 3 (5y((a + p){m + ~ k) + Me ^7 )] ' mdx ^'^ 

tl ~\~ fJ IsJUi LJJjo L/Jb% O L/Jjfc 

< CS\\(Vu,Vm)\\l 

- I u-Vm-mdx<C||u||3||Vm|| 2 ||m||6<C(5||Vm||^ (2.24) 

Jr 3 

! T -^Ve-Va-e&x<C\\Ve\\ 2 \\Va\\ 2 \\e\\ L ~<C5\\{Va,Ve)\\l, (2.25) 
J R 3 (a + py 

= — CjlX 

(a + p)(m + k) 

Cl£2 ^ [Me(|^ + ¥-) - k((« + P)(m + k) (2.26) 



(a + p)(m + fe) 9xj 9xj 9xi 3 
+ /*e^)]dx<C*||(V«,Ve)||l, 

-^ydx < C||£|| 2 ||e|| 6 ||e|| 6 ||m + fc|| 6 < C«J|| Ve||l, (2.27) 

in + k 



- [ u.Ve.eda;<C||u||3||Ve||2||e||6<C<J||Ve||l. 

JR 3 

The estimates (|2TT7|l - (|2~2g|l together with (|2TT4]) - (|2TT5)| imply 



(2.28) 



1 d 


\h\ 
m 

N 


2 + 

2 ' 
11 + 
11 + 


1 




2dV 

Id 


a + 
1 


P 


2dt" 
1 d 


a + 
1 


P 


2dt 


a + 


P 



|V/i||! <C5||(Vo,Vu)||i, (2.29) 

|Vm|||<C(J||(Vo,Vw)||l, (2.30) 

|Ve||i < C<5||(Va,Vw)|| 2 . (2.31) 

Step 2: L 2 -norms of V 3 u, V 3 /i, V 3 m, V 3 e 

Applying the differential operator di m n to the momentum equation of (|2.3j) . then 
multiplying it by di mn u, and integrating over R 3 , one gets 

2A l|ftmntt|1 ' 

1 /" . 1 

<9z m „( — — Aw)<9z mn ud:r + / 9j mn ( — — Vdivu)d; m „uda; 

a + P Jr 3 a + p , 

^ (2.32) 

dimn(u ■ Vu) ■ di mn udx + I di mn [f (a + p) - f(p)] ■ di mn divudx 

JR 3 

2 f 1 
7T / ^/m«[ — — V((a + p)(m + k))] ■ di mn udx. 

3 Jr 3 a + p 
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The first term on the right-hand side of ()2.32j) can be estimated as 



hmn( — —Au)di rnn udx 
a + p 

1 ,„ _ l2 . 



-\di mn Vu\ dx + - — — — -Vdi mn u : (Va <g> d lmn u)dx 
a + p J M3 [a + py 

6 2 

— -^diad m ad n aAu + - (di m ad n aAu + di n ad m aAu 



+ d mn adiaAu + diad m ad n Au + diad n ad m Au + d n ad m adiAu) [2.66) 

- t —K(di mn aAu + di m ad n Au + di n ad m Au + d mn adiAu 

[a + pY 

+ diad mn Au + d m adi n Au + d n ad hn Au)] ■ di mn udx 

< C5||(V 2 u,V 3 u,V 4 u,V 2 a,V 3 a)||^- / —L-\di mn V u\ 2 dx . 

Jr3 a + p 



Similarly, we can estimate the second and third term on the right-hand side of 
(j2"32) as 



dimn{ — —Vdivu)di mn udx 
a + p 



<C<5||(V 2 u,V 3 M,V 4 u,V 2 a,V 3 a)||^- / — !— \d lmn divu\ 2 dx, 

Js? a + P 



(2.34) 



di mn (u ■ V«) ■ d lmn udx < C8\\V 3 u\\l (2.35) 
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Now, let's estimate the fourth term on the right-hand side of (|2.32[) as 

dimn[f(a + p) - f(p)] ■ dimndivudx 
- / di mn [f(a + p)~ f(p)]di mn (— — — )dx 

[f'"{p)diad m ad n a + f"(p)d hn ad n a + f"(p)di n ad m a 

6 

(aTpY 



+ f"(p)d mn adia + f'{p)d lmn a] > - j- ■ diad m ad n a(a t + u ■ Va) 



=r^{di m ad n aa t + di n ad m aa t + d mn adiaa t + diad m ad n a t 

(a + pY 

diad n ad m a t + d n ad m adia t + u ■ V 'adi m ad n a + u ■ V 'adi n ad m a 

u ■ V ad mn adia + diu ■ V 'ad m ad n a + d m u ■ ~S/adiad n a + d n u ■ V 'adiad m a 

u ■ diX7ad m ad n a + u ■ <9 m Va<9/a9„a + u ■ d n V 'adiad m a) 
1 



-, 9 -(a t di mn a + dia t d mn a + d m a t di n a + d n a t di m a + di m a t d n a 
(a + pY 

+ di n a t d m a + d mn a t dia + u ■ Vadi„ ln a + diu ■ Vad mn a 

+ d m u ■ Vai9;„a + d n u ■ X7adi m a + u ■ V 'diad mn a + u ■ V<9 m a<9/„a 

+ u ■ Vd n adi m a + di m u ■ Vad n a + di n u ■ Va<9 m a + d mn u ■ Va<9;a 

+ u ■ V 'di m ad n a + u ■ V<9;„a<9 m a + u ■ Vd mn adia + diu ■ V9 m ai9„a 

+ diu ■ ~S/d n ad m a + d m u ■ ~S/diad n a + d m u ■ ~S/d n adia + d n u ■ ~S/diad m a 

+ d n u ■ Vd m adia) -\ —(di mn a t + di„ ln u ■ Va + di m u ■ V<9„a 

a + p 

+ di n u ■ Vd m a + d mn u ■ V<9;a + d L u ■ Vd mn a + d m u ■ V<9j„a 
+ d n u ■ Vd lm a + u ■ Vdi mn a)}dx 

<C<5||(Va,V 2 a,Va t ,V 3 a,V 2 a t ,V 2 u,V 3 u)||^- / [/'" (p)diad m ad n a 

Jr 3 

+ f" (p)di m ad„a + f"(p)di„ad m a + f" '(p)d mn adia + f'(p)di mn a) 
{di mn a t + u ■ Vdi mn a)]dx. 



(2.36) 



a + p 
Since 



1 /" 1 

f (p)dtad m ad n ad lmn a t dx = / <9;( — —f'"(p)diad m ad n a)d mn a t dx 



a + p Ji{3 a + p 

<C6\\(Va,V 2 a,V 2 a t )\\ 2 2 
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1 f 1 

f" (p)di m ad n adi mn a t dx = / di{ — —f (p)di m ad n a)d mn a t dx 



3 a + p J R 3 a + p 

<CT||(VyV 2 a t ,V 3 a)||2, 



1 f 1 ... ■ 

f"'(p)diad m ad n au-Vdi mn adx = / d j (— T ^f / "(p)diad m ad n au : >)di mn adx 



a + p Jm.3 a + p 



<C5||(V 2 a,V 3 a)||? 



and 



/ (p)d[ m ad n au-Vdi mn adx = / dj( — — / (p)di m ad n au>)di mn adx 



a + p J R 3 a + p 

<CS\\V 3 a\\ 2 2 , 



we obtain 



/ [f'"(p)diad m ad n a + f"{p)di m ad n a + f"{p)di n ad m a + f"(p)d mn adia] 

JH3 



x [ — —{di mn a t + u ■ Vdi mn a)} 
a + p 

< C5\{Va,V 2 a,V 2 a t ,V 3 a>" 2 



2- 

(2.37) 



Finally, 

f'JP) 
a + p 



di mn a(di mn a t + u ■ Vd lmn a)dx 



-If ^-.{d lmn a) 2 t dx+ l - f lifL(d lmn afdivu+^(d lmn a) 2 u-Va 

2 Jr3 a + p 2 J R 3 a + p a + p 

;(di mn a) 2 u ■ Vadx 



(a + f>) 2 

^-^ / ^T--(dimna) 2 dx+)- [ [H£l] t (d lmn a) 2 dx + C6\\V 3 a\\ 2 
2 dt J R 3 a + p 2 J R 3 a + p 

^-Hf —-{d lmn a) 2 dx + C5\\V 3 a\\l 
2 dt J M 3 a + p 



together with (|2.37[) . thus (|2.36[) can be replaced by 



dl mn [f{a + p) - /(/>)] • dimndivudx 

< C5\\{Va,Vat,V 2 a,V 2 a t ,V 3 a,V 2 u,V 3 u)\\ 2 2 (2.38) 

(di mn a) 2 dx. 



ld f /'(/"),« ^2 



2 dt J R 3 a + p 



/ 7n I l? 

I d mn [ — -^Va + Vm] • du mn udx 
Jr3 a + p 



GLOBAL SOLUTION FOR TURBULENT FLOW EQUATIONS 11 

The last term on the right-hand side of (|2.32|) can be estimated as 

2 f 1 _ - 

~o 5(m ™[ — —'V((a + p)(m + k))]-di mn udx 

3 Jk3 a + P 

2 f „ r m + k 
= 3 

— o / ["7 — , -^ d m ad n aS/a - - — — rr, ((m + k)d mn aVa + d m ad n mS7a 
3J R 3 (a + p) 3 (a + Pr 

+ d n ad m mVa + (m + /c)<9 m a<9„Va + (to + k)d n ad m Va) H = (<9 mn ?7jVa 

a + p 

+ d m mVd n a + d n mVd m a + (m + /c)Vd mn a) + <9 mn Vm] • du m nudx 

< C(5||(Va,V 2 a,V 2 77i,V 3 77i,V 4 m,V 4 u)||^ + C f (S7d mn a) 2 dx, 

Jm. 

which together with ([233]) - ([235)1 and (|235)) gives 

~[||ftmnU||2+ / tAPl(d hnn a) 2 dx}+ -^—\d imn V u\ 2 dx 

I dt J R3 a + p J R s a + p 

:\di mn divu\ 2 dx 



a + P (2.39) 

< Ctf||(Va, Va t ,V 2 a, V 2 a t , V 3 a, V 2 u, V 3 m, V 2 m, V 3 m, V 4 m)j|^ 

+ C / (V<9 m „a) 2 dx. 

Applying <9; mn to the energy equation, fc-equation and e-equation of (|2.3p . mul- 
tiplying the resulting equations by di mn h, di mn m and di mn e respectively, then 
integrating them over the whole space R 3 , one can prove that 
1 d in mi2 

2j t \\dlmnh\\ 2 

1 f *? 

dimni — —Ah) ■ dimnhdx + di mn ( — — ) ■ di mn hdx (2.40) 

> a + p J R3 a + p 

di mn (f (a + p)(a + p)divu) ■ dimnhdx - / di„ m (u ■ V/i) • di mn hdx, 
1 d 



2di 



|5j m „m||2 



= / d[ mn ( — —Am) -dimnindx- / di mn e ■ d lmn mdx 

JR3 a + p J R 3 

+ / dimni — — ) • di mn mdx - I di mn (u ■ Vm) • <9; mn md:r, 

j R 3 a + p J R3 



(2.41) 



1 <* „ 

Jlmn^- 1 1 2 



r-^ll<9; m .„e|| 2 



2di 



i re Cf 

dimn(-, — —-^Ae)-dim n edx+ di„ in (- --. =r) • di mn edx (2.42) 

{a + P) Jm 3 (a + p)(m + k) v > 



C F 2 f 

di mn { — tt) ■ dimnedx - / 5; m „(u ■ Ve) ■ di mn edx. 

R 3 m + k J R 3 
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As same as the estimate ()2.33j) . we can deduce 

dimn{ — —Ah)di mn hdx 

a + p 



<C<5||(V 2 /i,V 3 /i,V 4 /i,V 2 a,V 3 a)||^- / —^\di mn Vh\ 2 dx, 

Jr3 a + p 

f 

<C(5||(V 2 m,V 3 m,V 4 m,V 2 a,V 3 a)||^- / — i-^|<9 /m „Vm| 2 dir, 

Jk.3 a + p 



dimn( — —Am)di mn mdx 
a + p 



dimn( — —Ae)di mn edx 
a + p 

<CT||(V 2 e,V 3 £,V 4 £,V 2 a,V 3 a )|| 2 - / -^—\d lmn Ve\ 2 dx. 



(2.43) 



(2.44) 



(2.45) 



/R3 a + p 
Again, from the Holder inequality, (P^ ]) -(|2"^ ]) . W^ 2 (R 3 ) <-► L 6 (R 3 ) and H 2 (R 3 ) 

L°°(R 3 ), we get 



dlmn{ —) ■ dl mn hdx 

a + p 

d { l \( (— + —)--$■■ dyk ) dul Mt dp da 1 ] d hdx 

m ™ a + p dxj dxi 3 lJ dxk dxj (a + p) 2 dxj dxj m " 



2- 



< C6||(V 3 a, V 3 m,V 4 w, V 3 /i,V 4 /i)|| 

/ d imn ( — — ) • di mn mdx 

Js.3 a + p 

f du l du l dvP 2 - du 

= / di mn { — — - — [/z e (- \- ——)- -Sij({a + p)(m + k)+ p e - — )]} • di mn mdx 

J R 3 a + p oxj l oxj oxi 3 v axk 

< CS\\ (V 2 a, V 3 a, V 3 m, V 4 ?ti, V 3 u, V 4 u)|| 2 , 

9imn(-, =77 fz) ■ di mn edx 

R 3 (a + p)(m + k) 

_ r Ci£ <9w\ ,9m 1 Su- 7 , 2 , - 

,^^ (a + p)(m + fc) ^^ ( ^- + fe- ) "3^( (a + p)(m + fc) 

+ Me7j )]} ■ dlmnSdx 

OXk 

< C<5||(V 3 a,V 3 M,V 4 u,V 3 TO,V 3 £)|| 2 , 
d lmn (f'(a + p)(a + p)divu)-d lmn hdx < C6\\(V 3 a,V 3 u,V 4 u,V 3 h,V 4 h)\\l 

d lmn (u -Vh) ■ d lmn hdx < C6\\(V 3 u,V 3 h,V 4 h)\\ 2 , 

3 

/ di mn e ■ di mn mdx < C5\\ (V 3 £, V 3 m, V 4 77l) || 2 , 

■/K3 

^ m „(u-Vm)-9 im „mda;<C , 5||(V 3 m,V 4 TO,V 3 M)|| 2 , 
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dimni-^r) ■ d lmn edx < C<5||(W V 3 e)|| 2 , 
m + k 

- ! d lmn (u-We)-d lmn edx<C\\(W 3 e,W i e,W 3 u)\\l 

JR 3 

Incorporating the above estimates and ()2.40[) - (|2.42[) yields that 

\di m „h\\l+ / ——\di mn Vh\ 2 dx 

Jr 3 a + p 



Id , „„ f 1 

2~aV 

< C6\\ (V 2 a, V 3 a, V 2 h, V 3 /i, W 3 u, W 4 u)\\ 2 2 , 



(2.46) 



7;-T;\\ d lm n m\\2+ / — —\di„ m Vm\ 2 dx 

2 at J M3 a + p (2.47) 



< CS\\ (V 2 a, V 3 a, V 2 m, V 3 m, V 3 e, V 3 u, V V)| 



3„ V73„, Y74„,M|2 



(2.4 



(2.49) 



o'3lll 5 «mne||i+ / — —\d lmn \7e\ 2 dx 
2 at J R 3 a + p 

< CS\\ (V 2 a, V 3 a, V 3 m, V 2 e, V 3 £, V 3 u, V 4 w)|| \. 
Step 3: L 2 -norms of Va, V 3 a 
We first estimate for Va. For this purpose, we calculate as 

JR 3 * * 

= / Va ■ Va t + (a + p)a t Wa ■ u -\ Va t ■ u -\ Va ■ u t dx. 

Jr 3 2 2 

The first term of the right-hand side of (|2.49[) can be estimated as follows: 

/ Va ■ \7atdx = — / Va • Vdiv[(a + p)u]dx 

Jr 3 Jr 3 

= - / Va • (V 2 a • u) + Va • (Va • Vu) 

Jr 3 

+ (a + p)Va • Vdivu + |Va| divudx 

<C5\\Va\\\- / (a + p)Va • Vdivuda:. 

Jr 3 

Also, we estimate the rest three terms as 

(a + p)a t Va ■ udx < C5\\(a t ,Va)\\ 2 2 < CS\\(Va, Vu)\\l, (2.51) 

f (a + p) Va t -Mda;<CT||(Va,Vit)||^ (2.52) 

Jr 3 2 



(2.50) 
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^2 



(a + p) 

Va • u t dx 

2 

(a + p) 2 , 1 A 1 

Va ■ ( — u ■ Vm H -Au H = Vdivu 

2 a+p a+p 

- V[/(a + p) - /(/>)] - — J— V((a + p)(m + fc)))dx (2.53) 

<C5\\(Va,Vu)\\t- [ (a + ^ m + ^ |Va| 2 dx 

/■ ( a + P) 2 ^( p )| Va | 2 dx+ / ^Va • (Au + Wdivu)dx. 
Js.3 2 J R 3 2 



Since 



Va • (Am — Vdivw)da; = / ^a ■ (djjU 1 — didjV?)dx 

2 Jr3 2 

—diadjadjU 1 dijadjU 1 H — diadjadiU- 1 H dijadiU J 'dx = 0, 

together with (f!T4l?]) - (f2"33")) and ([23)1 - (j2l)|) . we show that 

/ [^|Va| 2 + ( a + d Va • u] t + C||Vo||l < C5\\Vu\\ 2 . (2.54) 

Jk 3 2 2 

Now, we turn to estimate for V 3 a. Almost parallel to the inequality (|2.49p . we 

get 

^\d hn na\ 2 + { «^d lmn ad lm u%dx 

-d lmn adi mn div[(a + p)u] + (a + p)a t di mn adi m u n (2.55) 

(a + p) 2 [a + p) 2 
-z di mn div[(a + p)u\dimU H di mn adi m u t dx. 

We estimate the right-hand side of (|2.55[) as follows: 
-di mn adi mn div{(a + p)u]dcc 

3 

/ -di mn a[u ■ Vdi mn a + diu ■ Vd mn a + d m u ■ V<9z„a + d n u ■ V5 im a 

JR3 



(2.56) 



+ d lm u ■ V<9„a + di„u ■ V<9 rn a + d mn u ■ \7d L a + di mn u ■ Va + di mn adivu 
+ di m ad n divu + di n ad m divu + d mn adidivu + diad mn divu + <9 m a<9z„divu 
+ <9 n a<9 im divu + (a + p)9; m „divu]da; 
< C<5||(V 3 a,V 3 M)||2 - / (a + p)di mn adi mn divudx, 

JR 3 

(a + p)a t d lmn ad lm u n dx < C6\\(V 3 a,X7 2 u)\\i (2.57) 
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2 

(a + p) 2 



di mn div[(a + p)u]di m u n dx 

di m u n [di mn {u ■ Vo) + dimnddivu + di m ad n divu 



+ di n ad m d\vu + d mn adidivu + diad mn divu + d m adi n drvu + d n adi m divu ^ ' ' 

+ (a + p)di mn divu]dx 

< CS\\ (V 2 a, V 3 a, V 2 m, V 3 u, V 4 u)\\j, 



2 Oi mn adi m u t dx 

( a + P)\ a r i a n 1 « „ 1 « < 
^ — T —^di mn adi rn [-u l d i u n H —d H u n H —d in u % 

I a + p a + p 

- d n (f(a + P)- f{p)) - ^—zd n ((a + p){m + k))]dx 

3(a + p) (2.59) 

< CS\\ (V 2 a, V 3 a, V 2 u, V 3 w, V 3 m, V 4 m)||^ 

di mn a{di m uu n + dimmU l )dx - / f'(p)(di mn a) 2 dx 

(a + p)(m + fc) 



3 

Noting that 



(di mn a) 2 dx. 



(a + /9)dj mn ad imn divuda; + / di mn a(d[ mil u n + di mm u l )dx 

JR3 * 

di mn a(di mii u n - di mni u l )dx 

•3 A 

-- diadi mn adi ml u n + (a + p)di mni ad hnl u n - d l adi mn adi mn u l 

1 JR3 

(a + p)di mm adi mn u l dx 
1 



diadi mn adi mi u n - d t adi mn adimnU l dx 
<Ctf||(V 3 a,V 3 «)|& 

based on the estimates (|2.55l) - (|2.59p . with the help of (|2.5|) - (|2.6[) and the interpola- 
tion inequality, one obtains 

[dd;„ m a| 2 + {a + p > di mn ad lm u n ] t dx + C||V 3 a||| 
* * (2.60) 

< C(5||(V 2 a, V 2 u, V 3 u, V 4 u, V 3 m, V 4 m)||^. 
Step 4: Conclusion 
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Consequently, multiplying (|2.39|) by a appropriate small constants a, together 
with ([235]) . ([2T2"g|) - (|2"3Tj) . (gg5J) - ([2Ilgl) , (gUS and (EUni), we have 



I™"' 



,V 3 /i,rn,V 3 m,£,V 3 e)||2+ / F(a)da; + a[||V 3 u||2 

JR3 



- / ^(d lmn a) 2 dx}+ f [i|Vo 



p (a + p) 2 

2 + -^Vfl'M 



! Ia ,2 (« + P)% n nlj ! ( 2 - 61 ) 

+ -|d im „a| H di mn adi m u \dx\ 

+ C(a)\\ (Va, V 3 a, Vu, V 4 w, V/i, V 4 /i, Vm, V 4 m, Vs, V 4 e)||l 

<0, 

where we have used the fact that 

i+i 
\V l a t \<Cj2(\V k a\ + \V k u\), i = 1,2. 
fc=i 

Integrating the inequality (|2.6ip . from (|2.4[) - (|2.6[) . (|2.12|) and the smallness of a, 
we can finish the proof of Proposition 12.11 □ 



3. Proof of global existence 



We will finish the proof of Theorem 1 1.1 1 in this section. First, let's state the local 
existence. Since it can be proved in a standard way as that in |27 | l29j. we omit the 
proof. 

Proposition 3.1. Under the assumption of the Theorem \1.1[ then there exists 

a constant T > such that the system \2.S\) admits a unique smooth solution 

(p, u, h, k, e) which satisfies that there exists a constant C% > 1 such that for any 

t€[0,T], 

r* 

2 i / /I|V7„I|2 i ll/Y7„. V7U Y7L T7^M|2 



\\(p-p,u,h,k-k,e)\\ 2 H3 + / (||Vp||^ + ||(V U ,V/i,Vfc,V £ )||^)cfe 

Jo (3.1) 

< C 2 \\(p - P,u ,h ,k - k,e )\\ 2 H3 . 

In the following, by a continued argument, combining the local existence and 
the a priori estimates proposition, we will prove the global existence of smooth 
solutions. 

First, suppose 



^o = ||(/5o - p,ua,h ,k Q - fc,£ )||fl-3 < mm(S/y/C^,6/ y / C 1 C 2 ), (3.2) 



GLOBAL SOLUTION FOR TURBULENT FLOW EQUATIONS 17 

where 5 is defined in Proposition 12.11 Since the initial data satisfy Eq < S/^/U^, 
then by Proposition 13.11 there exists a constant T* > such that there exists a 
unique solution on [0,T»] satisfying 

Ei := sup \\(p - p,u,h,k-k,e)\\ H 3 < ^/QzEq. (3.3) 

0<t<Ti 

Therefore, using the inequality Eq < d/\fCxC^, from Proposition 12. 1[ we have 

E x < y/C^E < 5/^C~ 2 . (3.4) 

Notice that T* depends only on Eq. Starting from T*, then the initial problem 
(|2.3[) with initial data (p, u, h, k, e)(T*) still has a unique solution on [T*,2T*], and 
from Proposition 13. 11 we get 



sup \\(p-p,u,h,k-k,e)\\ H3 < ^/C~ 2 E l < ^C l C 2 E < 6. 

T*<t<2T* 



Again from Proposition 12. 1\ one can deduce 

E 2 = sup \\(j> — p,u,h,k — k,e)\\ H 3<y/CiEo<$/y/C2- 

a<t<2T, 

Repeating the procedure for < t < NT*, N = 1,2,3,- ••, we can extend the 
local solution to infinity as far as the initial data are small enough such that Eq < 
min(<5/ v / C^, 5/ \fCiC 2 )- Thus the proof of Theorem 1 1.1 1 is complete. □ 
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